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Using data from observations made with XMMNewton, we present an X-ray 
analysis of two flares observed in £ Boo. The flare loop parameters are de- 
rived using various loop models including state-of-the-art hydrodynamic flare 
model. The loop lengths derived for the flaring loop structure are found to be 
less than the stellar radius. The exponential decay of the X-ray light curves, 
and time evolution of the plasma temperature and emission measure are sim- 
ilar to those observed in compact solar flares. The X-ray light curve of post 
^ ' flare phase is investigated with wavelet analysis. Wavelet analysis clearly show 

' oscillations of the period of 1019 s. Using the observationally estimated loop 

£\J , length, density and magnetic field, the theoretically derived oscillation period 

for fast-kink mode approximately matches with the observationally estimated 
■ period. This is the first likely observational evidence of fundamental fast-kink 

mode of magnetoacoustic waves in the stellar loops during the post-flare phase 
of heightened emission. 
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1. Introduction 



Sudden enhance and gradual decay in the intensity of a star is the most 
;_j ' known feature of a flare observed in it. During the flare large amount of 

energy is released in a short interval of time. Flares are known to be a man- 
ifestation of the reconnection of magnetic loops, accompanied by particle 
beams, chromospheric evaporation, rapid bulk flows or mass ejection, and 
heating of plasma confined in loops. Stellar flares could be radiating several 
orders of magnitude more energy than a solar flare. The flare produced by 
the stellar sources (e.g. RS CVn and BY Dra types) are usually detected 
in UV and X-rays. Flares from these stars present many analogies with 
the solar flares. However, they also show significant differences, such as the 
amount of energy released. Modeling the dynamic behavior of X-ray flares 
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allows us to constrain the properties of magnetic loops in ways that can- 
not be done from analysis of quiescent coronae that provide only a spatial 
and temporal average over some ensemble of structures. Therefore, study 
of stellar flares is a valuable tool for understanding stellar coronae, as these 
are dynamical events, which therefore embody different information than 
available from static, time-averaged observations (e.g., 1 , 2 , ) . 

With the advent of HINODE / STERIO the study of oscillations in the 
solar corona became the most interesting now a days. The waves in the solar 
corona have been explained by magneto-hydrodynamic (MHD) oscillations 
in coronal loops. 5 MHD waves and oscillations in the solar and probably in 
the stellar atmospheres are assumed to be generated by coupling of complex 
magnetic field and plasma. These MHD waves and oscillations are one of 
the important candidates for heating the solar/stellar coronae and accel- 
erating the supersonic winds. In the Sun, magnetically structured flaring 
loops, anchored into the photosphere, exhibit various kinds of MHD oscil- 
lations (e.g. fast sausage, kink and slow acoustic oscillations). Stellar MHD 
seismology is also a new developing tool to deduce the physical properties 
of the atmosphere of magnetically active stars, and it is based on the anal- 
ogy of the solar coronal seismology. 6 One of the most exciting aspects of 
"coronal seismology" is that it potentially provides us with the capability 
for determining the magnetic field strength in the corona, 7 as well as, in the 
stellar case, with information on otherwise unresolved spatial scales, e.g., 
flare loop lengths. 

In this article, we have studied the X-ray flares and associated MHD 
oscillations in the post flare phase of the star £ Boo. £ Boo is a nearby (6.0 
pc) visual binary, comprising a primary G8 dwarf and a secondary K4 dwarf 
with an orbital period of 151 yr. In terms of the outer atmospheric emission, 
UV and X-ray observations show that the primary dominates entirely over 
the secondary. 8 Recently, using the Chandra observations Wood & Linsky 9 
show that £ Boo primary contributes 88.5% of total X-ray emission. 

2. Observations and data reduction 

£ Boo was observed with XMM-Newton satellite using the EPIC MOS 
detectors and Reflection Grating Spectrometer (RGS) on 2001 January 19 
at 19:25:06 UT for 59 ks. The details on the XMM-Newton satellite and 
its detectors are given in the references 10 and. 11 The data were reduced 
using the Science Analysis System version 8.0 with updated calibration 
files. Details of data reduction is given in. 2 The stellar separation of xi Boo 
primary and xi Boo secondary is 6.3 arcsec, which more than the spatial 
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resolution of XMM-Newton. Therefore, we can safely consider that £ Boo 
primary is observed by XMM-Newton. 



3. X-ray light curves 

The MOS light curve of £ Boo in the energy band 0.3-10 keV is shown in 
Figure 1. Two flare like features, marked as Fl and F2 are clearly seen. 
The quiescent state in the light curve is marked by Q. A higher emission 
marked by U is also seen in the X-ray light curve. Both flares are found to 
be ~ 2.6 times more energetic than that of the quiescent state. A flux of 
2 x 10 -11 erg cm~ 2 s _1 was observed during the quiescent state of £ Boo. 
The e-folding decay times of the flares Fl an F2 were found to be 4396 ±333 
and 9885 ± 362 s, respectively. However, the e-folding rise times were found 
to be 3169±259 and 3351 ±244 s for the flares Fl and F2, respectively. The 
duration of flares Fl and F2 were 8.4 and 11.3 ks. From the light curve, 
it is seen that before ending the flare Fl the flare F2 starts. Similar loop 
systems have been observed in a flare on the Sun (e.g. so-called Bastille Day 
flare 12 ), and in a stellar analogue dMe star Proxima Centauri. 4 In these two 
events, a double ignition in nearby loops was observed or suggested, and 
the delay between the ignitions appears to scale with the loop sizes. Both 
flares are appeared to be a long decay flare (jd > 1 li). 



Time (s since JD 2451928.99394) 

Fig. 1. X-ray light curve of the star £ Boo detected by MOS1, vertical lines are different 
time segments. 
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4. Evolutions of flares 

In order to trace the spectral changes during the flares, we have analysed 
the spectra of the different time intervals shown in Fig. 1. Flare Fl and 
F2 have binned into five and six intervals, respectively. The spectra corre- 
sponding to different time segments associated with flares Fl and F2 are 
shown in Fig. 2 (a) and 2 (b), respectively. The quiescent state spectra 
are also shown by solid circles in each Fig. 2 (a) and 2(b) for comparison. 
The quiescent state spectra of £ Boo were fitted with a single (IT) and 
two (2T) temperature collisional plasma model known as APEC, 13 with 
variable elemental abundances. Only 2T plasma model with variable abun- 
dance (Z < 0.23Zq) were found acceptable. The temperature and emission 
measure (EM) were found to be 0.20 ±0.05 keV and 15 ±6 x 10 51 cm~ 3 for 
cool component, and 0.57 keV and 12 ± 1 x 10 51 cm -3 , respectively. The 
luminosity during the quiescent state was found to be 5 x 10 28 erg s _1 . 




0.5 1 2 5 0.5 1 2 5 

Energy (keV) Energy (keV) 

(a) (b) 



Fig. 2. (a) Spectra of different time segments during the flare Fl, and (b) Spectra of 
different time segments during the flare F2. The quiescent state spectra is also shown by 
solid circles in each Figure (a) and (b) 

To study the flare emission only, we have performed IT spectral fits of 
the data, with the quiescent emission taken into account by including its 
best-fitting 2T model as a frozen background contribution. This is equiva- 
lent to consider the flare emission subtracted of the quiescent level, allows 
us to derive one effective temperature and one EM of the flaring plasma. 
The abundances were kept fix to that of the quiescent emission. Figs. 3 (a) 
and 3 (b) show the temporal evolution of the temperature and correspond- 
ing EM of the flares. Both the temperature and EM show the well-defined 
trends i.e. the changes in the temperature and the EM are correlated with 
the variations observed in the light curves during the flares. The hardness 
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ratio can reveal the information about temperature variations. The hard- 
ness ratio also varied in the similar fashion to their light curves (see also. 2 




Fig. 3. (a) Evolution of temperature, and (b) emission measure during the flare Fl and 
F2 and (c) the densitytemperature diagram, where EM 1 / 2 has been used as a proxy of 
density. 



5. Determinations of density using RGS spectra 

Figs. 4 (a), 4 (b) and Fig. 4 (c) show the He-like triplets from O vn during 
the quiescent state, flare state Fl and flare state F2+U. Line fluxes and po- 
sitions were measured using the XSPEC package by fitting simultaneously 
the RGS1 and RGS2 spectra with a sum of narrow Gaussian emission lines 
convolved with the response matrices of the RGS instruments. The contin- 
uum emission was described using Bremsstrahlung models at the tempera- 
tures of the plasma components inferred from the analysis of EPIC-MOS1 
data. The emission measure derived from the analysis of the EPIC data 
was used to freeze the continuum normalization. The detailed study of the 
RGS spectra was done by Nordon & Behar. 14 Therefore, for the present 
purpose, we give the line fluxes of He-like triplets from O vn in Table 1. 




(a) (b) (c) 

Fig. 4. RGS spectra of £ Boo during quiescent state Q (a), flare state Fl (b), and flare 
state F2 (c). 
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Tabic 1. Observed O VII line fluxes during the Q, Fl and F2+U 
states 

Quiescent (Q) Flare (Fl) Flare (F2+U) 

A(A) Flux a A(A) Flux" A(A) Flux" 



22.08(f) 2.1 ±0.2 22.07 2.3 ± 0.2 22.08(f) 2.2 ± 0.3 
21.78(i) 0.8 ±0.2 21.77 0.9 ± 0.2 21.77(i) 0.9 ± 0.2 
21.58(r) 3.0 ±0.3 21.56 3.5 ± 0.3 21.57(r) 3.4 ± 0.3 

Note: a Measured flux in 10~ 4 photons/cm 2 /s 

The He-like transitions, consisting of the resonance line (r) Is 2 x So 
- Is2p 1 pi, the intercombination line (i) Is 2 1 So - ls.2p 3 P1, and forbid- 
den line (f) Is 2 1 So - Is. 2s 3 Si are density- and temperature-dependent. 
The intensity ratio G = (i+f)/r varies with temperature and the ratio R 
=f/i varies with electron density due to collisional coupling between the 
metastable 2 3 S upper level of forbidden line and the 2 3 P upper level of 
intercombination line. In the RGS wavelength, the O vn lines are clean, 
resolved and potentially suited to diagnose electron density and tempera- 
ture (see Fig. 4). These lines at 21.60A(r), 21.80A(i) and 22.10A(f) have 
been used to obtain temperature and density values from the G- and R-ratio 
using CHIANTI database (version 5.2.1; 15 ). The G-ratio as function of tem- 
perature and intensity ratio R as a function of density are shown in Figs. 
5(a) and 5(b), respectively. It is clear from Fig. 5(a) the temperature was 
found to be ~ 2 MK. Therefore, for the representative temperature around 
2MK, we derive electron densities of n e = 1.3+q' 2 x 10 10 cm -3 for quiescent 
state (Q) and 1.0±J;| x 10 10 cm" 3 during flare Fl and 1.6±£;§ x 10 10 cm" 3 
for flare state (F2+U). The electron densities were found well with in a er 
level for each three states. These densities are found to be similar to that 
of solar conditions. 

6. Loop modeling 

In a star, flares cannot be resolved spatially. However, by an analogy with 
solar flares and using flare loop models, it is possible to infer the physical 
size and morphology of the loop structures involved in a stellar flare. Based 
on quasi static radiative and conductive cooling during the early phase of 
the decay. Haisch 16 suggested an approach to model a loop. Given an esti- 
mate of two measured quantities, the emission measure (EM) and the decay 
timescale of the flare (jd ), this approach leads to the following expression 
for the loop length {L} la )\ 
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Fig. 5. (a) The ratio G = (i + f )/r of the summed intensities of the OVo intercombi- 
nation and forbidden lines over the intensity of the recombination line (A= 21.60 A) and 
(b) intensity ratio R = f /i of the OVu forbidden (A = 22.12 A) and intercombination 
(A = 21.80 A) lines as a function of electron density, calculated using the CHIANTI 
database. 



L ha (cm) = 5 x W-^EM) 1 ^ 4 (1) 

After Haisch 16 various stellar flare models came up. These are (i) two- 
ribbon flare method, 17 (ii) quasi-static cooling method, 18 (iii) pure radia- 
tion cooling method, 19 (iv) rise and decay time method, 20 and (v) hydrody- 
namic method, 2122 The two-ribbon flare model assumes that the flare decay 
is entirely driven by heating released by magnetic reconnection of higher 
and higher loops and neglects completely the effect of plasma cooling. The 
other four methods are based on the cooling of plasma confined in a single 
flaring loop. The hydrodynamic model includes both plasma cooling and 
the effect of heating during flare decay. Reale et al. 21 presented a method 
to infer the geometrical size and other relevant physical parameters of the 
flaring loops, based on the decay time and on evolution of temperature and 
the EM during the flare decay. The thermodynamic loop decay time can be 
expressed 23 as 



3.7 x I(T 4 L 
T th = ^ (2) 

V max 

where L is the loop half-length in cm, and T max is the flare maximum 
temperature (K) calibrated for EPIC instruments as, 

Tmax = 0.13T ofw . 6 (3) 
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Tobs is the maximum best-fit temperature derived from single temperature 
fitting of the data. The ratio of the observed exponential light curve decay 
time Td to the thermodynamic decay time r t h can be written as a function 
which depends on the slope £ of the decay in the density-temperature plane. 
For the XMMNewton EPIC spectral response the ration Td/r t h — F(() is 

99 

given as 

— = F ^ = 7~~ttW + L36 (for 0.35 < C < 1.6) (4) 

T th C - 0.35 

Combining above equations the expression for semi-loop length is 



L ~ Lhyd ~ 3.7 x 10- 4 F(C) 1 j 

Alternatively, Reale 22 derive the semiloop length from the rise phase 
and peak phase of the flare as 

T 5/2 

L hyr = 3xW 5 /H M ^ (6) 

Tm and tu are temperature and time at which density peaks. 

Based on magneto-hydrodynamic simulations Shibata & Yokoyama 24 
assumed that, to maintain stable flare loops, the gas pressure of the evap- 
orated plasma must be smaller than the magnetic pressure 

2nkT < g (7) 

where B is the minimum magnetic field necessary to confine the flaring 
plasma in the loop. The derived equations for B and L p are 



„/ EM \ 3/5 / n c \-V5/ T V 8/5 



The average temperatures of the loops, usually lower than the real loop 
maximum temperatures, are found from the spectral analysis of the data. 
According to equation (3), the loop maximum temperatures for the flare 
Fl and F2 are found to be 18.9 ±0.8 and 25 ±2 MK, respectively. The peak 
values of temperature are similar to those found in the solar flares. However, 
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Tabic 2. Loop parameters obtained from various loop models. 



Flare 


Lha 


Lhyd 


Lhyr 


Lp 


B 


n e 




(10 10 cm) 


(10 10 cm) 


(10 10 cm) 


(10 10 cm) 


(G) 


(10 1(, cm- 3 ) 


Fl 


2.2 ±0.2 


1.9 ±0.8 


1.7 ±0.2 


5.7 ±1.4 


18.3 ±3.4 


1 0+ 1 ' 5 


F2 


3.8 ±0.2 


7.9 ±0.6 


3.6 ±0.8 


3.6 ±0.9 


29.4 ±6.0 


1 6+ 1 - 2 
1 -°-0.8 



the peak emission measure was found to be an order more than those found 
in the solar flares. The path of density-temperature (n-T) diagram is shown 
in Fig. 3(c). The EM 1 / 2 was taken as a proxy of density. The solid lines 
represent the best linear fit to the corresponding data, providing the slope 
£. The values of £ for the flare Fl and F2 were found to be 0.7 ± 0.3 and 
2.1 ±0.4, respectively. This indicates that the flare Fl is driven by the time- 
scale of the heating process, whereas sustained heating is negligible during 
the decay of the flare F2. 

The semi- loop lengths derived from above methods are given in table 2. 
The loop length derived for flare Fl is consistent using three methods (see 
equations (1), (5) and (6)). However, the loop length derived from pressure 
balance method ( see equation 9) was found to be ~ 3 times more than 
other three methods. The value of £ for the flare F2 is outside the domain 
of the validity of the hydrodynamic method (see equation5), therefore, the 
loop length derived from hydrodynamic method based on decay phase may 
not be actual one. The loop length derived from other three methods (see 
equations (1), (6) and (9)) is found to be consistent at a value of 3.6 x 10 10 
cm. The magnetic field derived using the equation (8) is 18.9 and 29.4 Gauss 
for the flares Fl and F2, respectively. 

The total energy rate is defined as H = Eh x V, where Eh = 
10~ 6 T^ X L~ 2 erg s _1 cm -325 heating rate per unit volume and V = 
EM/n 2 volume of the loop. Therefore, the total energy rate was deter- 
mined to be 5.4 x 10 30 and 9.8 x 10 30 ergs s _1 for the flare Fl and F2, 
respectively. The peak X-ray luminosity during the flare Fl and F2 was 
found to be 1.02 x 10 29 and 1.11 x 10 29 ergs s _1 , which is only 2 and 11 % 
of the total energy rate, respectively. 

7. Search for coronal oscillations 

We have used the wavelet analysis IDL code "Randomlet" developed by 
E. O'Shea. The program executes a randomization test, 2627 which is an 
additional feature along with the standard wavelet analysis code 28 to ex- 
amine the existence of statistically significant real periodicities in the time 
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series data. The advantage of using this test is that it is distribution free 
or nonparametric, i.e., it is not constrained by any specific noise model 
such as Poisson, Gaussian, etc. Using this technique, many important re- 
sults have been published by analyzing approximately evenly sampled data 
(e.g.,, 27 , 2930 ). The wavelet power transforms of the U part of the X-ray light 
curve of cadence 40 s are shown in the left panel of Fig. 6, where the darkest 
regions show the most enhanced oscillatory power in the intensity wavelet 
spectrum. The crosshatched areas are the cone of influence (COI), the re- 
gion of the power spectrum where edge effects, due to the finite lengths 
of the time series, are likely to dominate. The maximum allowed period 
from COI, where the edge effect is more effective, is 7363 s. In our wavelet 
analysis, we only consider the power peaks and corresponding real periods 
(i.e., probability > 95%) below this threshold. The period with maximum 
power detectable outside the COI is 2883 s with a probability of 94% and 
the repetition of only 3 cycles in the time series. However, another peak is 
visible at a period of 1019 s with a probability of 99% and the repetition 
of 11 cycles. Therefore, only the periodicity of 1019 s satisfies the recently 
reported stricter criteria of O'Shea & Doyle, 31 i.e., at least four cycles of 
repetition over the lifetime of the oscillations. The right panel in Fig. 6 
shows the global wavelet power spectra of the time series from which the 
statistically significant period (1019 s) is selected. Wavelet analysis was also 
performed in the decay phase of the flare F2, and no significant periodicity 
was found. 




50 100 150 200 250 300 0.001 0.010 0.100 

Time (minutes) Power 



Fig. 6. The wavelet power spectrum of the U part of X-ray light curve of £ Boo. 
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8. Evidence of Fast-kink Oscillation 

In this section, we probe the presence of various MHD modes associated 
with the flare F2 observed in £ Boo. The overdense magnetic loops are pres- 
sure balanced structures and may contain fast kink and sausage oscillation 
modes with phase speed greater than the local Alfvenic speed (w/k > va), 
and slow acoustic modes with sub- Alfvenic speed (co < va) under coronal 
conditions. The expressions for the oscillation periods of slow (P s ), fast- 
kink (Pfk), and fast-sausage (Pf s ) fundamental modes are given as follows 

( 32 733-j . 

2L 2L 
P s = — = — 

JC T J C 

Pi ^ M ^j'\ iA ^ s (12) 

where L, a, p and B are loop length, loop width, ion mass density and 
magnetic field, respectively. The ct — c va/{c 2 + v A ) X ^ 2 , where c Q and va 
are sound and Alfvenic speeds respectively. The Ck is phase speed. Under 
coronal conditions, B Q w B e , p e « p a , nn ~ n e , and c Q << va- The 
subscripts ' d and 'e' stands for 'inside', and 'outside' of the loop. The j is 
the number of nodes in fast-kink and slow oscillations, and is equal to 1 for 
the fundamental mode. The approximate relation in equations (10) - (12) 
are given in coronal conditions for the fundamental oscillation periods of 
different modes in magnetic loops. Therefore, these approximate equations 
are also valid for the stellar loops in the coronae of magnetically active 
Sun-like stars. The parameters are expressed in units of L\q = L/10 1Q cm, 
a 8 = a/10 8 cm, n 9 = n o /10 9 c m - 3 , T = T e /10 e K, and B w = B/IQG. 

Using observationally estimated loop length (3.6 ± 0.8 x 10 10 cm) 
and width(a = 0.1L), density (1.6^8 x ^ W cm_3 ) an d magnetic field 
(29.4 ± 6.0G) for flare F2, and equations (10)-(12), the estimated oscil- 
lations periods P s , Pfk and Pf s for fundamental modes are found to be 
1586 ± 353, 1004 ± 391 and 313 ± 121 s, respectively. The observationally 



Lin 
1300^-s 

Vt 6 



(10) 
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derived oscillation period ~ 1019 s approximately matches with the theo- 
retically derived oscillation period for fundamental fast kink mode. These 
oscillations may be formed either by the superposition of oppositely propa- 
gating fast-kink waves or due to flare generated disturbances near the loop 
apex. The flaring activity may be the possible mechanism of the generation 
of such oscillations in the stellar loops. Cooper et al. 34 have found that the 
transverse kink modes of fast magnetoacoustic waves may modulate the 
emission observable with imaging telescopes despite their incompressible 
nature. This is possible only when the axis of an oscillating loop is not 
orientated perpendicular to the observer's line of sight. Using this theory 
O'Shea et al. 35 have also found the first two harmonics in kink oscillations 
in cool solar loops. Similar effect may also cause the modulation in the emis- 
sion from stellar loops. Resonant absorption is the most efficient mechanism 
theorized for the damping of a kink mode in which energy of the mode is 
transferred to the localized Alfevanic oscillations of the inhomogeneous lay- 
ers at the loop boundary. 36 Intensity oscillations due to the fast-kink mode 
of magnetoacoustic waves do not show any decay of the amplitude during 
the initial 250 minutes in the post-flare phase observations of £ Boo. This is 
probably due to an impulsive driver, which forces such kind of oscillations 
and dominates over the dissipative processes during the flare F2 of £ Boo. 

9. Summary 

The observed flares in G-dwarf £ Boo are similar to the solar arcade flares, 
which are as strong as those in M-dwarfs and are much smaller than the flare 
observed in dMe star, giants and pre-main-sequence analogous. The loop 
lengths are found to be smaller than the stellar radius. The electron densi- 
ties during the flares and quiescent state were found to be of the order of 
10 10 cm -3 . We also report the first observational evidence of fundamental 
fast-kink oscillations in stellar loops during a post-flare phase of height- 
ened emission in the cool active star £ Boo. Based on the solar analogy, 
such observations may shed new light on the dynamics of the magnetically 
structured stellar coronae and its local plasma conditions. 
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